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We study a class of one-dimensional interacting particle systems 
^v^j , with random boundaries as a microscopic model for Stefan's melting 

and freezing problem. We prove that under diffusive rescaling these 
particle systems exhibit a hydrodynamic behavior described by the 
■ solution of a Cauchy-Stefan problem. 

I 1. Introduction. In this work we return to the classical Stefan's freezing 

on the ground model [8]. It could be described in the following way: Consider 
the real line occupied by a heat-conducting material (heat is transmitted 
only by conduction). This material is initially almost everywhere character- 
. ized by a bounded and measurable temperature function T : K — > R. Accord- 

J> | ing to the temperature the material could be in one of two phases, a liquid 

phase for positive temperatures and a solid phase for negative temperatures. 
The temperature T = is that of crystallization at which both phases may 
If} ' occur. The problem consists in determining the temporal evolution of the 

temperature profile. 

We consider this problem under more restrictive conditions. Suppose that 
at initial time the liquid phase fills the domain u > at positive tempera- 
tures and the solid phase fills the domain u < at negative temperatures. 
Denote by p ^ : M_ — > M_ and p\ : M + — > M + the initial temperature profile. 
We are able to determine a function B = B(t) describing the time evolution 
of the boundary between the two phases and their temperature functions, 
respectively and p\{t,u) for the solid and liquid phases. It is well 

known that these functions satisfy a Cauchy-Stefan problem: 



dtp-i = a^duup-t, d t pi = aid uu pi, 
(1.1) B(t) = k{a„ 1 d u p„ 1 (t,B(t))-a 1 d uPl (t,B(t))}, B(0) = 0, 
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Pi (t,B(t)) = 0, Pi (0,-)=ft U (-), 

where p°_i : K_ — ► K_ and p1:R+ are bounded measurable functions, 

a_i > and ai > are the coefficients of heat conduction of the material 
with respect to the solid and liquid phases and k > is a scaling factor for 
the temperature. 

In this paper we present a microscopic model for Stefan's equation through 
an appropriated interacting particle system and scaling limit techniques. 
Such sorts of descriptions have been proposed previously by Chayes and 
Swindle [2] in the case of finite domains with a_i = 0. Rezakhanlou [6] and 
Bertini, Butta and Riidiger [1] also derive Stefan's equation as hydrodynamic 
limit of interacting particle systems. 

We shall denote by Z, N and Z_ the sets of integers, positive integers and 
nonpositive integers respectively, and by K_ and M+ the sets of nonpositive 
and nonnegative real numbers. 

For the informal description of the microscopic model, consider the one- 
dimensional lattice Z with each site being occupied by a molecular agglom- 
erate of type —1 for the material in the solid state and of type 1 for the 
material in the liquid state. According to its internal energy, each agglomer- 
ate is classified by a heat unit of or 1. An interaction between neighboring 
sites occurs independently in the following way: If the particles are of the 
same type, then their heat units are interchanged after a mean a_i expo- 
nential time for particles type —1 and after a mean a\ exponential time for 
particles of type 1. If the particles are of distinct type and their heat units 
are also distinct, at rate 1 the heat unit of the agglomerate whose heat value 
was 1 drops to and simultaneously the other agglomerate changes type. If 
the particles are of distinct type and their heat units are equal to 1, both 
heat units drop to after a mean 1 exponential time. Moreover, we start 
with configurations such that the agglomerates are of type —1 if they are at 
the left of the origin; otherwise they are of type 1. 

We will show that this system has a hydrodynamic behavior under diffu- 
sive scaling described by the solution of a Cauchy-Stefan problem of type 
(1.1) with scaling factor k = 1, where the temperature is the macroscopic 
heat density profile. The general case with an arbitrary k can be obtained 
from the previous one rescaling the temperature by k~ l . Here the diffusive 
scaling is expected since the hydrodynamic behavior of the simple symmet- 
ric exclusion process is described in this scale. Actually, our model could 
be described as a coupling between two one-dimensional nearest-neighbor 
simple symmetric exclusion processes in the semi-infinite lattice. To make 
this identification, consider each agglomerate of solid (resp. liquid) phase as 
a site in the space of the Z_ (resp. N) in such a way that this association 
preserves the order. At each site whose associated agglomerate has heat unit 
equal to 1 we put a particle. In each of the lattices Z_, N particles evolve as 



STEFAN'S MELTING AND FREEZING PROBLEM 



3 



in a nearest-neighbor one-dimensional simple symmetric exclusion process 
with jump rates a_i and a±, respectively. Superposed to this dynamics, a 
particle at the boundary of one of the lattices waits a mean 1 exponential 
time and attempts to leave the system. If no particle occupies the bound- 
ary of the other lattice, this particle vanishes triggering a translation of the 
whole system to the right or left depending whether the particle was occu- 
pying a site in Z_ or in N. If the boundary of the other lattice is occupied, 
then both particles leave the system. 

From the technical point of view, the main difficulty of this problem lies 
in the fact that no entropy argument can be used due to the annihilation 
mechanism. Indeed, if one fixes the boundary to be at the origin, the unique 
invariant measure is degenerate, and even when estimating the relative en- 
tropy with respect to a nonstationary state as in [3, 5], the translations 
introduce expressions too large to be estimated by the sole entropy. Cou- 
pling is therefore the unique available tool in this context. 

2. Main result. Let T = {(—1,0), (—1, 1), (1,0), (1, 1)} Z be the configura- 
tion space and denote a typical configuration by (a, rj) = {(a(x), r](x))} x< =z G 
T. Fix a_i > 0, a\ > and define the generators Go-, Gi by 

@o/)M) = E iM*) = < x + i)W(*){f(<rrf' x+1 ) - fMh 

ri) = £ l{a(x) + a(x + l)}{f(T^ +1 (a, r,)) - f(a, r,)}. 

xez 

In these formulas, / stands for a cylinder function f :T — ► K, rj x ' y for the 
configuration 7] with spins at x, y interchanged, 

!rj(x + l), if z = x, 
V (x), i£z = x + l, 

7](z), otherwise, 

and T x ' x+1 (a, rj) = {<J,r)) for the configuration defined as follows: 

!cr(x), if z = x + 1, rj(x + 1) = 0, rj(x) = 1, 

a(x + 1), if z = x, r)(x) = 0, rj(x + 1) = 1, 
cr(z), otherwise, 

_/ s_ JO, if z = x,x + 1, 

\rj(z), otherwise. 

Hence, two neighboring particles of different type annihilate each other at 
rate 1: (—1, 1), (1, 1) — > (—1, 0), (1, 0); and a particle sitting next to a vacant 
site of different type dies transforming the neighboring site into a site of its 
type at rate 1: (-1,0), (1,1) -> (1,0), (1,0). 

Denote by (at,rjt) the Markov process associated to the generator Q = 
Go + Q\ speeded up by iV 2 . The goal of this article is to show that its 
macroscopic behavior is described by solutions of Stefan's equations. 
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The Cauchy-Stefan problem. Let a(u) = a_il{n < 0} + a\\{u > 0}, 
A(u) = ua(u). Fix a bounded measurable function pq : R — ► R and consider 
the Stefan problem 

d t p = a{p)Ap, 

(2.2) ui(t) = bi{d u A(p(t,Ui(t)-)) - d u A(p(t )Ui (t)+))}, 

p(o,0 = po(0- 

Here A stands for the Laplacian, U{{t) for the curves at which p(t,Ui(t)) = 
0, h = l{p(t, Ui (t)-) < < p(t, Ui (t)+)} - l{p(t, Ui (t)+) < < p(t, Ui (t)-)} 
and the first equation should be understood throughout [0, T] x R, except 
on the curves Ui(t). 

Denote by Cq' 2 {[0,T) x R) the set of functions G:[0,T)xR-tl with 
compact support which are continuously differentiable with respect to the 
first variable and twice continuously differentiable with respect to the second 
variable. A bounded measurable function p : [0, T] x R is said to be a weak 
solution of the Stefan problem (2.2) if for every function G G Cq' 2 ([0, T) x R), 

-oo rT 

du j dt {A(p(t, u))AG(i, u) + u(p(t, u))d t G(t, u)} 

(2.3) 



oo JO 



/+oo 
duuj(p (u))G(0,u) = 0, 
-oo 

where 

P<0, 

(2.4) uj(p) = {p, p>0, 




p = 0. 

The proof of uniqueness of weak solutions for the Stefan equation (2.2) 
presented in [7], Theorem 20, page 312, for boundary- valued problems can 
be easily adapted to our context. Furthermore, the generalized solution is 
continuous according to Theorem 21 of [7]. 



a{z) 



The initial states. Let A be the subset of T of all configurations (<r, rf) 
for which there exists x in Z such that 

— 1, if z < x, 
1, if z > x. 

Note that A is stable under the dynamics induced by the generator Q. 
For every (a, if) G A, let 

b = b(a) := sup{z : a(z) = — 1} 

be the boundary of the configuration (<r, 77). 

In the proof of the hydrodynamic behavior of the process (cr, 77) we impose 
some conditions on the initial states. Let {m N ,N > 1} be a sequence of 
probability measures on T. We assume that: 
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(HI) For every N > 1, m is concentrated on configurations of the set A 
such that b(a) = 0. 

(H2) There exists a bounded measurable function po : K — ► [—1, 1] such that 



Po(«) ^ > 0, 



Po(u) du < 



for all a > and such that for each continuous function G : 
compact support and each 5 > 0, 



lim to 

TV— >oo 



TV 



5^ G(x/JV)ff(x)T7(a 



xG2 



duG(u)po(u) 



>5 



with 



0. 



Notice that condition (HI) forces the initial profile to be negative on R_ 
and positive on M + . 

In the case a_i = 0, we impose one more condition: 



(H3) For every N>1, 



m N {(a, r/)eA: n(x) = 0, x < 0} = 1. 



JTie hydrodynamic behavior. For each probability measure m on T con- 
centrated on A, denote by the probability measure on the path space 
.D(R+,r) induced by the Markov process (a t ,r] t ) with generator Q speeded 
up by N 2 and initial measure m. 

Denote by M. = A4(K) the space of signed Radon measures on R en- 
dowed with the vague topology. Integration of a function G with respect to 
a measure it in M. is denoted by (tt, G). To each configuration (rj,a) £ F 
we associate the empirical measure 71"^ = ir N (r],a) in M by assigning mass 
a(x)N" 1 to each particle: 

7lJV = Tf^2 a ( x )v(x)5 x /N- 
Letn?=ir N (a t , Vt ), bf = b(a t )/N. 

Theorem 2.1. Fix a sequence of initial measures {m N : N > 1} satisfy- 
ing assumptions (HI), (H2) and (H3) i/a_i = 0. For eac/i i > 0, as N f oo, 
i/ie empirical measure ir^ converges in probability to an absolutely continu- 
ous measure n(t,du) = p(t,u)du, whose density p(t,u) is the weak solution 
of the Stefan problem (2.2): For every continuous function G:R — >R with 
compact support and every 5 > 



lim I 

TV^oo 



<7rf,G>- / duG(u)p(t,u) 



>S 



0. 
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Moreover, for every 5 > 

lim F f m N N [\b(t)-B(t)\>5} = 0, 

N— >oo 

where B is the solution of B(0) = 0, 

(2.5) B(t) = a. 1 (d uP )(t,B(t)-)-a 1 (d u p)(t,B(t)+). 

3. The fixed boundary model. Recall the definition of the set A given 
in the previous section and of the boundary b = b(a, rj) of a configuration 
in A. To each configuration (a, rj) in A let £ = in = {0, 1} Z be the 
configuration viewed from the boundary: 

= + 6). 

It is not difficult to check that £t is a Markov process with generator C given 
by 

C = a-\L\ + Lb + aiL 2 . 

Here L\ and L2 are the parts of the generator related to the motion of 
particles in a simple symmetric exclusion process on Z_ and N respectively: 

L\ = L x -i^ x , L2 = ^ Lx,x+ii 

x<0 x>l 

where, for every local function / : £1 — > R and every integer x, 

(L x , x+1 f)(o=f(e> x+1 )-f(o, 

and ^ :c ' a; + 1 is the configuration £ with spins at x, y interchanged defined in 
(2.1). 

In contrast, Lf, is the part of the generator related to the dissipative 
feature of the system: For every local function / : — ► R 

(£&/)(£) = £(1)[1 - e(0)]{/(r_!(£ - ft)) - /(£)} 

+ ^(o)[i-e(i)]{/(ri(£-eo))-/(0} 

+£(o)e(i){/(e-^o- qi) -my, 

where g x stands for the configuration with no particles but one at x, and 
{t x :i£Z} for the group of translation so that {t x ^){z) = £(z + x) for all z 
in Z. 

Fix a bounded measurable function Ao:R— >R such that \§(u) > 0, 
Ao(^) < for n > 0, n < 0, respectively. A pair (A, -D), where A is bounded 
measurable function A:[0,T]xl^M strictly positive a.e. on (0, 00), strictly 
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negative a.e. on (—00, 0) and D : [0, T] — > R is a bounded variation continuous 
function vanishing at t = 0, is said to be a weak solution of 

d t X = aAX + D(t)d u X, 

(3.1) D(t)=a. 1 (d u X)(t,0-)-a 1 (d u X)(t,0+), 
A(0,-) = A (-), 

in the layer [0, T]xK + , if for every function G G C^QO, T) x R), 

+00 /-T 

/ dt {A(\(t, u))AG(t, u) + A(t, u)d t G(t, u)} 
-00 Jo 

/•T H-00 

(3.2) ~ dt {X(t,u)d u G(t,u)-G{t,0)}dD(t) 

JO J-00 

/+00 
dtiA (n)G(0,n) = 0. 
-00 

Notice that we are requiring the solution to be strictly positive, negative on 
(0,oo), (—oo,0), respectively. 

Lemma 3.1 below shows that weak solutions of (2.2) may be obtained 
by a simple change of variables from weak solutions of (3.1). In particular, 
uniqueness of weak solutions of (3.1) follows from the uniqueness for (2.2). 
Indeed, assume that (Xt,Dt) is a weak solution of (3.1). By Lemma 3.1 and 
by the uniqueness of weak solutions of (2.2), p(t,u) = X(t,u — D t ) is the 
unique weak solution of (2.2). Since the weak solution of (2.2) is continuous 
and since A is a.e. strictly positive, negative on (0, 00), (— 00, 0), respectively, 
for each t > 0, p(t,-) vanishes at a unique point. This determines uniquely 
Dt and therefore A. 

Lemma 3.1. Let (Xt,Dt) be a weak solution of (3.1) and let p(t,u) = 
X(t,u — Dt). Then, p is a weak solution of (2.2). 

Proof. Consider a weak solution (Xt,Dt) of (3.1) and write Dt as the 
difference of two continuous increasing bounded functions: Dt = Df — D^ . 
Let D ±,e be smooth uniform approximations of D^ and set Df = D^' £ — 
D t ' e so that 

(3.3) lim sup \D £ t - D t \ = 0. 

£ ^°0<t<T 

Fix a smooth function G : [0, T] xl-*M with compact support and van- 
ishing at the boundary t = T. Let H £ (t,u) = G(t,u + Df). H £ is a smooth 
function with compact support. Therefore, since (Xt,Dt) is a weak solution 
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of (3.1) and since G vanishes at the boundary t = T, 
0= (X T ,G(T,u + D £ T )) = (X T ,H^) 

(3.4) = {X 0l H £ ) + I da (A s , (d 8 + aA)Hl) 

Jo 

- [ T {(X s ,d u H s s )-HI(0)}dD s . 
Jo 

Recall the definition of the function u> given in (2.4) and that weak solutions 
of (3.1) are strictly positive in (0,oo) and strictly negative in (— oo,0). The 
first two terms on the right-hand side of (3.4) can be rewritten as 

(lu(X ),H £ )+ [ ds(X s ,aAH £ s )+ [ ds (w(A a ), d s H e 3 ) 
Jo Jo 

+ (l{(-oo,0)},# E ) + [ T ds(l{(-oc,0)},d s H £ ). 
Jo 

Since G, and therefore H £ , vanish at the boundary t = T, the second line of 
the previous expression is equal to 0. On the other hand, the first two terms 
of the first line are easily seen to converge to 

(3.5) <w(po),G )+ [ T ds{A( Ps ),AG s ) 

Jo 

as e I 0. The last term of the first line together with the last term on the 
second line of (3.4) is equal to 

(w(A a ), (d s G)(s,u + D £ s )) ds + f T (u;(X s ), (d u G)(s,u + D £ s )) dD £ s 



f T {(X s ,(d u H £ )(s,u)) - H £ s (0)}dD s . 
Jo 



10 

The first term converges, as e [ 0, to 

(3.6) [ T ds{uj(p s ),d s G s ), 

Jo 

while, by definition of u and H £ , the sum of the second and third terms is 
equal to 

-{ G(s,D £ s )d(D £ -D s )+ [ (X s ,d u H £ )d(D £ s -D s ). 
Jo Jo 

It is not difficult to show from (3.3) that this expression vanishes as e | 0. 
It follows from (3.4), (3.5) and (3.6) that 

rp p 

(u( Po ),G )+ [ ds(A( Ps ),AG s )+ [ ds(u;(p s ),d s G s ) = 0, 



which concludes the proof of the lemma. □ 
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Hypotheses on the initial measures. Fix a sequence of probability mea- 
sures {fi N :N > 1} on Q. To prove the hydrodynamic behavior of the system 
we will assume that 



(HI) There exists a bounded measurable initial profile Ao 
that 



-1,1] such 



Xq(u) du > 0, 



Xq(u) du < 



for all a > and such that for each S > and each continuous function 
G : R — ► R with compact support 



lim 

N-^oo 



jV 



J2G(x/N)l ± (x)^ 



+oo 



du G(u)\q(u) 



>5 



0, 



where l±(u) = -l{u < 0} + l{u > 0}. 
In the case a_i = 0, we impose one more condition: 
(H2) For every N > 1, 

/{^fi:^)=0,i;<0} = l. 

T/ie hydrodynamic behavior. For each probability measure /i on !], de- 
note by P^ the probability measure on the path space Z)(R + ,0) induced 
by the Markov process £j with generator C speeded up by N 2 starting from 
the initial measure \x. 

Let D+(t) [resp. D^(t)] be the total number of particles on N (resp. 
Z_) which left the system before time t divided by N and let D N (t) = 
D N (t) - (t). Formally D+ (t) = N' 1 E*>i{£o(z) ~ 

Theorem 3.2. Fix a sequence of initial measures {[j, N ,N > 1} satisfying 
(HI) and (H2) if a_i = 0. Then, for any t > 0, any continuous G:R — ► R 
with compact support and any 5 > 



lim P^v 



l^l±(x/iV)G(x/iV)6(x)- | d«G(«)A(i, 



>6 



0. 



lim F» N [\D N (t)-D(t)\>5}=0, 



where (X,D) is the unique weak solution of (3.1). 



4. Proof of Theorems 2.1 and 3.2. We first show how Theorem 2.1 can 
be recovered from Theorem 3.2. 
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Proof of Theorem 2.1. Notice first that the evolution of the original 
process ((Tt,r)t) can be derived from the one of £t since = D N (t) and 

vt = T--6((T t )6- 

By Theorem 3.2, for every t > 0, converges in probability to D(t). 
This proves the second statement of the theorem since -D(-) satisfies (2.5) in 
virtue of (3.1). 

On the other hand, if G : R — > R is a continuous function with compact 
support, by the previous relations between £ and (r],a), 

= W E l±(x/N)G(D»(t) +x/NMx). 

By Theorem 3.2, this expression converges in probability to / duG(D(t) + 
u)X(t,u) and this integral is equal to / duG(u)p(t,u) by Lemma 3.1. This 
concludes the proof of Theorem 2.1. □ 

We turn now to the proof of Theorem 3.2. We present the proof in the 
case a_i > which is more difficult. The same arguments apply to the case 
o_i = 0. 

Recall that we denote by M = A4 (R) the space of signed Radon measures 
on R endowed with the vague topology and that we denote integration of 
a function G with respect to a measure tt in M by (tt,G). For each iV > 1 
and each configuration £ of $7, let 7r^ be the empirical measure associate to 
£ given by 

Note the indicator function l±(a;) in the definition which corresponds to 
considering particles on Z_ as having negative charge. Let tt^ = Tt N {£,t) and 
recall that we are speeding up the process by N 2 . 

Recall that D+(t) [resp. D^(t)] stands for the total number of particles 
on N (resp. Z_) which left the system before time t divided by ./V and that 
D N (t)=D^(t)-Df(t). 

With this notation, Theorem 3.2 states that the sequence Q^iv converges 
weakly to the probability measure concentrated on paths (irt,D(t)) whose 
first coordinate is absolutely continuous n(t, du) = X(t, u) du, the density 
being the solution of (3.1) (cf. [4]). The proof consists in showing tightness, 
that all limit points are concentrated on absolutely continuous paths which 
are weak solutions of (3.1) and uniqueness of weak solutions of this equation. 

Uniqueness of weak solutions of (3.1) was discussed in the previous section, 
while tightness is proved at the end of this section. We show now that all 
limit points are concentrated on weak solutions. 
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Fix a sequence // of probability measures on f2 satisfying the assumptions 
of the theorem. Note that all limit points of the sequence = Q^jv are 
concentrated on absolutely continuous measures since in the limit the ir- 
measure of a finite interval is bounded by its Lebesgue measure. 



Proposition 4.1. All limit points Q* of the sequence Q^ N are concen- 
trated on trajectories (irt,Dt) such that 

- (ir ,G) = f t ds(7r s ,(d s + aA)G s )- f{{v s ,d u G s ) - G(s,0)}dD s 
Jo Jo 

for every t>0 and G in Cq' 2 {[0,T) x R). 



Proof. The proof of this proposition is divided in several steps. We 
start examining some martingales associated to the empirical measure. Fix 
G £ Cq' 2 ([0,T) x R) and 5 > 0. Consider the martingale M G ' N given by 



Ml 



G.N 



(Ttf ,G t ) - {4, G ) - f\d s + N 2 £)(^,G s )ds. 
Jo 



An elementary computation shows that the quadratic variation (M G,N )t of 
this martingale is equal to the time integral of 

1 



w J2 a(x/NMx + 1) - t(x)] 2 {(V N G)(x/N)Y 



(4.1) 



x^O 

+ £(0)[1 - m{G(-l/N) + (tt n , (V N G)(x - l/N))} 2 
+ am-m){G(2/N) - (n N ,V N G)} 2 
+ e(0)£(l){G(l/A0-G(0)} 2 . 



In particular, since G is a smooth function with compact support, by Cheby- 
shev and Doob inequalities, 



sup \M t ' | > 

0<t<T 



<45~ 2 E" N [{M C 



< C(a, G)5~ 2 1 ^ + Efr y\cs (0) + 6(1)} ds 



for some finite constant C(a,G) depending only on a(-) and G. Therefore, 
by Lemma 5.5, 



(4.2) 



lim P^v 



sup \M^ N \ > 5 

0<t<T 



0. 
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On the other hand, an elementary computation shows that for every 
smooth function H :R — > R with compact support, 

N 2 C(7r N ,H) = (Tr N ,aA N H) 

+ N[£(l) - ,V N H) - N£(1)H(2/N) 

(4.3) + N£(0)H(-1/N) + [oie(l) + a_if (0)](VjvfO(0) 

+ £(0)(7r JV ,A JV iT) 

-e(o)e(i){(7r JV ,A Jv iy)-iv- 1 (A A riy)(i/iv)}, 

where and Vat denote respectively the discrete Laplacian and gradient. 
Therefore, in view of Lemma 5.5, up to negligible terms, the martingale 
M t ,N can be written as 



(Trf , G) - (ttq ,G) - [ ds (Trf , d s G s + aAG s 
Jo 



dsN[Ul) -£ s (0)]{« AG S ) - G(s,0)}. 



By Lemma 4.3, provided we let e [ after N f oo, we may replace -/V[£ s (0) 
6(1)] bye- 1 [5 Ar (s + e)-D Ar (s)]. Therefore, in view of (4.2), 



lim lim sup 

£^0 n^oo 



sup 

0<t<T 



<7rf,G)-«,G> 
t 



_v 



/ ds(Tr?,(d s + aA)G s 
Jo 



t D N (s + e)-D N (s) 
~\~ I ds 



x d u G s )-G(s,0)} 
Hence, for any limit point Q* of the sequence Q , 

(7T t , G) - (tto, G) - f ds (tvs, (d s + aA)G s ) 
Jo 

D(s + e) - D(s) 



>6 



lim ( 

e-*0 



sup 

0<t<T 



+ 



ds ■ 



■{(7r s ,d u G s )-G(s,0)} 



>5 



0. 



By the proof of tightness of the second marginal of Q N presented at the 
end of this section, Q* is concentrated on paths D which are of bounded 
variation. In particular, if d e (s) = e~ 1 {D(s + e) — D(s)}, for any continuous 
function H : [0,T] — s-R, /q ds d £ (s)H(s) converges, as e { 0, to J H(s) dD(s). 
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Since by Lemma 4.6 Q* is concentrated on paths 7r< which are continuous 
for the vague topology, the proposition is proved. □ 

The second main result of this section states that limit points of the 
sequence Q N are concentrated on trajectories 7Tj whose density is bounded 
below by a strictly positive function. 

Proposition 4.2. For each 5 > 0, there exists a strictly positive con- 
tinuous function R$ :R — > (0, 1] with the following property. All limit points 
Q* of the sequence Q N are concentrated on trajectories (irt,Dt) such that 



n{i)> j R s {u)du 



for all <t <T and all finite intervals I = [c, d] such that c > 5. A similar 
statement holds in (— oo,0). 

The proof of this proposition is postponed to Section 5. 

Lemma 4.3. Fix a smooth function G in Cq' 2 ([0, T] x R) and 5 > 0. 
Then, 



lim lim sup" 
£-+0 tv->oo 



j>7V 







sup 


fA 


_0<t<T 





D%(s + e)-D%(s 



(vrf ,G S ) 



>5 



0. 



A similar statement holds if we replace £ s (l), D+(s) by £ s (0), D^(s), re- 
spectively, or (tt^,G s ) byG(s,0). 

Proof. A simple computation shows that 



M»(t) = D»(t)-N Ul)ds 
Jo 

is a martingale with quadratic variation (M+)t given by Jq £ s (1) ds. We may 
therefore write e~ 1 {D I ^(s + e) — D+(s)} as 



Mf (s + e)-Mf (s) 



+ N 



S+£ 



&(l)dr. 



The martingale part is easy to estimate because it vanishes in the limit 
N | c<3. Indeed, by Chebyshev and Schwarz inequalities and by the explicit 
formula for the quadratic variation of the martingale (t), 

M^{s + e)-M^(s) 



IsJV 







sup 


f ds 


,0<t<T 


Jo 



-(ir^,G s ) 



C(G) N 



T+e 



ds\M%(s)\ 



< C(G 2 T)_ N 



>S 

r rT+e 







1/2 
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Here and below C(G), C(G,T) are finite constants depending only on G 
and G, T, respectively. By Lemma 5.5, this expression vanishes as N f oo. 

It remains to consider the difference iV£ s (l) — N J^ +e £ r (l) dr, which is 
slightly more demanding. We first perform a time integration by parts to 
obtain that 



^ds|iVe s (l)-iv£ +E ^(l)dr|{7rf,G. 
<C(G)|iV jf d«&(l) + N J* +£ 



+ 



dsNUl)\(^,G s )-- I 



Lemma 5.6 permits to estimate the first term on the right-hand side. To 
estimate the second one, write the difference (tt^,G s } — (Tr^,G r ) as 



M?> N - M^ N + f\d v + N 2 £)(ir^,G v ) dv. 

Jr 



On the one hand, 



(4.4) 



sup 

0<t<T 



dsNUV 1 



sup I M t 

0<t<T+e 



e 



G.N \ 



dr{M, 



G,N 



s—e 
T 



>8 



o 



N£,(l)d8>6/2 



Fix 7 > 0. By Lemma 5.5, there exists a finite constant A for which 

rT 

' <7 



/ N£.(l)d8>A 
Jo 



for all N > 1. Therefore, the previous probability is less than or equal to 



7 + 



sup \M?' N \>6/2A 

0<t<T+e 



It follows from (4.2) that this expression is bounded by 7 as N j 00. This 
proves that (4.4) vanishes in the limit N f 00. 
It remains to show that 



(4.5) 



sup 

0<t<T 



*dsNZ a {l)^ 



dr I (d v + N 2 C){rr^,G v )dv 



>5 



vanishes as N f 00, e j 0. 

By the explicit expression for (d s + N 2 £)(tt s , G s ) given in (4.3), we have 
that the absolute value of the integral in this formula is dominated by 

C(a,G) [ T NWL)dsle+ sup N[&0) +6(1)] da 

Jo I o<t<rJt 



STEFAN'S MELTING AND FREEZING PROBLEM 



15 



Repeating the argument presented just after (4.4) to eliminate J Q N£ s (l) ds 
and applying Lemma 5.6 to estimate the second term, we show that (4.5) 
vanishes as N f oo, e J, 0. This concludes the proof of the lemma. □ 

We conclude this section by proving that the sequence of probability mea- 
sures <Q> N is tight, which in our context reduces to showing that the marginal 
of on each coordinate is tight. We start with the empirical measure. De- 
note by Qf the marginal of on the first coordinate. 

Recall that is tight if for each smooth function with compact support 
G : R — > R, (tt^ ,G) is tight as a random sequence on J5(R+,R). Now fix such 
a function. To prove tightness for (nf* , G) it is enough to verify the following 
two conditions: 

(i) The finite-dimensional distributions of (tt^ ,G) are tight. 



(ii) For every 5 > 

lim lim sup ' 

£->0 tv->oo 



sup \{ir t ,G) - (ir a ,G)\ >5 

\s-t\<e 



0. 



Condition (i) is a trivial consequence of the fact that the empirical mea- 
sure has finite total mass on any compact interval. In order to prove condition 
(ii), consider the martingale with respect to T given by 

M?' N = (n?,G) - {ttq , G) - f N 2 C(^, G) ds. 

Jo 

Here the index N indicates that we are considering the process speeded up 
by N 2 . Therefore, 

(*? , G) - (vrf , G) = M t G > N - M?' N + jf* iV 2 £(vrf, G) dr. 

From the previous expression, condition (ii) is a consequence of the next two 
lemmas. 



Lemma 4.4. For every 5 > and every function G in C 2 
pact support, 

lim lim sup ¥^ N sup 

e^O tv^oo M [\ s -t\<e 



f N 2 £(ir?,G) dr 



>S 



with com- 



Proof. In view of (4.3), since G is a smooth function with compact 
support, the expression inside the absolute value is bounded above by 

H 



C(a,G)L + J N{£ r (0)+Z r (l)}dr 



for some finite constant which depends on o, G only. To conclude the proof, 
it remains to recall Lemma 5.6. □ 
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Lemma 4.5. For every function G in C 2 (M) with compact support and 
every 5 > 0, 



lim lim supP^jv 



sup \M®' N -. 

s-t\<e 



G,N\ 
s 



>8 



0. 



Proof. Denote by (M G ' N ) t the quadratic variation of the martingale 



G N 

M t ' . By the Doob inequality, 



sup \M t G ' N -M^ N \ 

\s-t\<e 



>8 



< 



sup \M G ' N \ >5/2 



0<t<T 



<^A(M g ' n )t}. 



By the explicit expression for (M G ' N }t given in (4.1), the previous expres- 
sion is bounded by 



T 



{6(0) + &(!)} <fe 



5 2 {N ' ~r 

To conclude the proof of the lemma, it remains to apply Lemma 5.5. □ 

We turn now to the tightness of the second marginal of Q N . Since D^(0) 
-D+ (0) = 0, we need only to show that 



lim sup lim sup P^v 



sup \D»(t)-D%(8)\>8 

S-t\<£ 







for every 5 > and a similar statement for D^(t) in place of D+(t). In fact, 
we claim that for every 5 > 0, there exists e > such that 



(4.6) 



limsupP^v 

N->oc * 



sup \D%(t)-D%(s)\>8 

\s-t\<e 



0. 



Since D+(t) is increasing, the previous probability is bounded above by 



Te~ 



J2V» N iD»(\j + l}e)-D»(je)>8/2]. 
i=i 

It follows from Lemma 5.4 of the next section that for each 5 > 0, there 
exists e > such that 



limsupP^W [£>£ (t + e)- D 1 ^ (t) >6]=0 

uniformly in < t < T. This proves that the second marginal of is tight. 

We summarize in the next lemma what we just obtained. Notice that we 
proved tightness in the uniform topology. 
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Lemma 4.6. The sequence Q is tight in the uniform topology. In par- 
ticular, all limit points are concentrated on continuous trajectories for the 
vague topology. 

Proof of Theorem 3.2. By Lemma 4.6, the sequence is tight in the 
uniform topology and by Propositions 4.1 and 4.2, all limit points are con- 
centrated on weak solutions (A, D) of (3.1). By uniqueness of weak solutions, 
presented in Section 3, the theorem is proved. □ 

5. Coupling. We prove in this section some important but technical re- 
sults which are used in the previous sections. The proofs rely on a coupling 
between the process Ct defined in Section 3 and an exclusion process Ct sim- 
ilar to Ci with the difference that the configuration is not translated when 
a single particle dies at the boundary. The generator £ of this process is 
therefore a_iL_i + a\L\ + L' b , where L' b is given by 

(4/)(C) = C(i)[i-C(0)]{/(C-ft)-/(C)} 
+ C(o)[i-c(i)]{/(C-f?o)-/(C)} 
+ C(0)C(i){/(C- 00-01) -/(C)}- 

Notice that both marginal processes on Z_ and on N behave as an ex- 
clusion process with disappearance at the boundary, whose hydrodynamic 
behavior is well known. The leading idea of this section is to show, through 
appropriate couplings, that the original process does not differ much in sev- 
eral aspects from the one defined above. 

Denote by {Q ■ t > 0} the Markov process with generator £ speeded up 
by N 2 and recall that we denote by {Ct : t > 0} the Markov process with 
generator C speeded up by N 2 . Let D^([s,t\) be the total number of £- 

particles in N which died in the time interval [s,t]. D7([s,t]), D^([s,t]) are 
defined analogously. 

Lemma 5.1. There exists a coupling (Ct,Ct) f or which 
(5.1) ( [0, t] ) + D+ ([0, t]) < 2D£ ([0, t]) + 2D+ ( [0, t] ) 

for all t > 0. 

Proof. The coupling between £t and Q can be described as follows. 
Assume that the initial configurations are identical at time 0: Co = Co- Label 
all particles and denote by Xf (resp. Y t 3 ) the position at time t of the jth £- 
(resp. C-) particle. We assume that X 3 Q < X§ if j < k, X§ < < Xq, X j q = Y(j 
for all j. 

The £- and C-partides with the same label jump together preserving the 
order of the labels until a particle dies. If two C-particles die simultaneously, 



18 



C. LANDIM AND G. VALLE 



we couple the disappearance of the £- and (-particles. If it is a single £- 
particle which disappears, assume, without loss of generality, that it has a 
positive label and denote this time by To. Note that due to the translation, 
X?p o = Yj, Q + 1 for all labels j associated to alive particles. 

Denote by T\ the first time after To in which the total number of disap- 
pearances of ^-particles in N is equal to the total number of disappearances 
of ^-particles in Z_: 

7i = inf{t > T : D+([T , t]) = D^([T ,t})}. 

In the time interval po,Ti], the coupling forces the £- and ("-particles 
with the same label to jump together. It may happen, however, that a (- 
particle on N disappears while its corresponding ^-particle remains alive. In 
this case, the ^-particle becomes a second-class particle to allow the coupled 
particles to jump together. The same phenomenon may occur on Z_, where 
a ^-particle may disappear while its corresponding (-particle remains alive. 
In this case also the (-particle becomes a second-class particle. Notice, in 
particular, that the difference Xf — does not depend on j for coupled 
particles. 

Due to the translations, for any Tq < t < T±, the total number of £- 
particles which died on N is bounded by the total number of (-particles 
which died on N: 

D+([T ,t])<D+([T ,t}). 

Also due to our definition of T\, for any To < t < T\, the total number of 
^-particles which died on Z_ is bounded by the total number of ^-particles 
which died on N: 

D^([T ,t])<D+([T ,t]). 

Therefore, (5.1) holds for any < t < T±. 

To proceed by iteration, let K\ = D^([Tq, T±]) be the total number of 
(-particles which died in N in the time interval po,Ti], and let k\ = 
D^([Tq,Ti\) be the total number of ^-particles which died in N in the time 
interval po,Ti]. By definition of the coupling, at time T±, there are K\ — k\ 
uncoupled ^-particles on N and K\ — k\ (-particles which died on N associ- 
ated to the Ki — k\ uncoupled ^-particles. These K\ — k\ (-particles should 
not be forgotten, since they will be used to compensate the eventual death 
of the second-class uncoupled ^-particles. The important fact for the recur- 
rence argument is that the number of dead (-particles which were not used 
to compensate the death of ^-particles is at least equal to the number of 
second-class uncoupled ^-particles. Notice also that there might be second- 
class uncoupled (-particles on Z_. Since they do not play any role in the 
argument, we do not refer to them again. 
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Assume that the first single ^-particle to die after T\ is in Z_. If it were in 
N, we could repeat the arguments presented in the last paragraph to arrive 
at the same conclusions obtained there and iterate again the argument. 

Denote by T 2 the first time after T\ in which the total number of disap- 
pearances of ^-particles in N is equal to the total number of disappearances 
in Z_: 

T 2 = inf { t > Tx : D+ ( [Tx , t] ) = DJ ([Tx , t] )}. 

The coupling in [Ti,T2] is the same described before, in which coupled par- 
ticles jump together until one of them dies. For Tx < t < T 2 , let L 2 (t) = 
D^([Tx,t}), L 2 = L 2 (T 2 ), £ 2 (t) = D^([Tx,t}), £ 2 =£ 2 (T 2 ). By definition of 
the coupling D^([Tx,t]) < l 2 {t) < L 2 (t), so that (5.1) holds for < t < T 2 . 
On the other hand, at time T 2 , there are: 

(a) at most Kx — kx uncoupled ^-particles on N. There might be less since 
a second-class ^-particle might have died, its death being compensated by 
the death of a (-particle on Z_, 

(b) L 2 — i 2 uncoupled ^-particles on Z_, 

(c) Kx — kx C-particles which died on N (in the time interval [To,Ti]) and 
which are associated to the remaining uncoupled ^-particles on N, 

(d) L 2 — t 2 ("-particles which died on Z_ associated to the remaining 
uncoupled ^-particles on Z_. 

Thus, at time T 2 , the total number of uncoupled second-class ^-particles 
is still smaller than the total number of dead and disassociated (-particles. 

Assume, without loss of generality, that the first single particle to die 
after T 2 is in N. Denote by T3 the first time after T 2 in which the total 
number of disappearances of ^-particles on N is equal to the total number 
of disappearances in Z_: 

T 3 = inf{i > T 2 : D+([T 2 ,t}) = DJ([T 2 ,t])}. 

The coupling remains the same. The next argument, though elementary, 
requires much notation. For T 2 < t < T3, let: 

( a ) J3(t) be the total number of second-class ^-particles which die on N 
in the time interval p2,£], 

(b) ks(t) — j3(t) be the total number of first-class ^-particles which die 
on N in the time interval [T 2 , t] , 

(c) Ks(t) — js(t) be the total number of (-particles which die on N in the 
time interval [T 2 , t] , 

(d) £s(t) be the total number of ^-particles which die on Z_ in the time 
interval [?2,t]. 
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By the definition of the coupling, £%(t) < ^(t) < K$(t). The js(t) second- 
class ^-particles which died on N were associated to (-particles which died 
before. Since there is a factor 2 in (5.1), we also associate to these (-particles, 
j3(t) A^3(t) ^-particles which died on Z_. The ks(t) — js(t) first-class £- 
particles which died on N are taken care of by ks(t) — J3(t) (-particles which 
died on N. The factor 2 in (5.1) allows to include (£ 3 (t) - h(t)) + < k 3 (t) - 
J3(t) ^-particles which died on Z_. Up to this point we showed that all 
disappearances of ^-particles in p2,i] can be compensated by disappearance 
of (-particles in [T2,t] and by disassociated ("-particles which died before T2. 
Therefore, (5.1) holds in the time interval [22^3]. 

To be able to iterate this argument notice that there are K^T^) — £3(13) 
second-class ^-particles created on N in the time interval [T2, T3] and 1^3 (T3) — 
^3(^3) (-particles which died in this interval and whose deaths were not used 
to compensate £ deaths. We may therefore associate these new second-class 
^-particles to these newly dead (-particles and iterate the argument. This 
concludes the proof of the lemma. □ 



Denote by (j - , ( 4 + the marginals of the process (j on Z_, N, respectively. 
Notice that both marginals evolve as an exclusion process in which particles 
leave the system at the boundary. This system plays an important role in the 
sequel and deserves a notation. For b > 0, denote by (3t the Markov process 
on {0, 1} Z+ with generator £ = £& given by 

(£/)(/?) = bY,{W x ' x+1 ) - /(/?)} + /9(0){/G9 - 00) - /(/?)}• 

For T > and a measure p on {0, 1} Z +, denote by P^ the probability 
on D([0, T], {0, 1} Z + ) induced by the Markov process 0t speeded up by iV 2 
starting from p. Expectation with respect to P^ is denoted by Ejj. 

It is well known that the process fit has a hydrodynamic description. Let 
tt^ be the empirical measure associated to (3t : vf^ = iV _1 Y^x>o Pt(x)5 x /N- 

Proposition 5.2. Consider a sequence of probability measures p N on 
{0, 1} Z + such that 



lim u 

iV^oo 



N 



(tt n ,G) - I Po (u)G(u)du 




>5 







for some measurable function po : M + — ► [0, 1] , every 5 > and every contin- 
uous function with compact support G. Then, 



lim ] 

N^oo 



(Trf ,G) - / p(t,u)G(u)du 




><5 
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for every 5 > and continuous function with compact support G, where p is 
the solution of the linear equation 

dtp = bAp on R+, 

(5.2) 

p(t,0)=0, p (0,-)=Po(-)- 

The proof of this result is similar to the one of Theorem 2.1 in [5]. More- 
over, the solution of (5.2) can be represented in terms of a standard Brow- 
nian motion Wt with absorption at the boundary u = 0: 

(5.3) p(t,u) = E u [p (VbW t )]. 

The coupling presented in Lemma 5.1 together with the hydrodynamic 
behavior stated in Proposition 5.2 permit to estimate the total number of 
particles which left the system in the original process This is the content 
of the next two lemmas. Recall the definition of D+(t), D^(t) introduced 
in Section 4. Denote by 1 the probability measure on {0, 1} Z + concentrated 
on the configuration with all sites occupied and by Dp(t) the total number 
of /3-particles which left the system before time t. 

Lemma 5.3. There exists a finite constant Co depending only on a\, a_i 
such that 

N^oo 

for all t>0 and all probability measures p. The lemma remains in force if 
is replaced by . 

Proof. By Lemma 5.1, the expectation in the statement of the lemma 
is bounded above by 2E^ [D^(t)], where D^(t) stands for the total num- 
ber of particles which left the system in the time interval [0, t] for the £ 
process. Since both marginals of £ evolve as an exclusion process with dis- 
appearance at the boundary, the previous expectation is bounded above 
by 4maxf )=aii0 _ 1 [Dp (i)]. By monotonicity, this latter expectation is less 

than or equal to 4max& =0lj0 _ 1 E±[Dp (t)]. By the hydrodynamic limit of /3, 
this expectation converges, as N j oo, to 

POO 

4 max / {1 — pb(t, u)} du, 
b=ai,a~i Jo 

where pt is the solution of (5.2) with initial condition po constant equal to 
1. With this initial condition, the solution of this equation can be written 
as pb(t,u) = 1 — 2P[Bt > u/V2b], where Bt is a standard Brownian motion. 
In particular, the previous displayed equation is equal to 

4V2t max{ y/a~{, ^ }E [| B x \ } . 

This concludes the proof of the lemma. □ 
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Lemma 5.4. For every 5 > 0, T > 0, there exists e > such that 
lim sup P^ [£>+(* + e) - D%(t) >5}=0 

uniformly for < t < T. The statement remains in force if we replace 
byD N . " 

Proof. Fix < t < T. Denote by fi N (t) the state of the process at time 
t. With this notation the probability appearing in the statement can be 
written as 

By Lemma 5.1 and by attractiveness of the /3-process, the previous expres- 
sion is less than or equal to 

2 max f?[D%(e)>5/4). 

6=ai,a_i 

By the hydrodynamic limit of f3 and the proof of the previous lemma, (e) 
converges in probability to 

fOO 



poo 

/ {1- p b (e,u)}du = V2beE[\B 1 
Jo 



In particular, if e is chosen small enough for the last expression to be less 
than 5/4, the previous probability vanishes as N ] oo. This concludes the 
proof of the lemma. □ 



Lemma 5.5. For every t>0, 



sup E at 

N P 



, t 

N 







W{6(0) + &(!)} da 



< oo. 



Proof. Recall the definition of the martingale M+{t) introduced in 
the beginning of the proof of Lemma 4.3. In particular, E^v [Jq N£ s (1) ds] = 
E^v [D^(t)]. It remains to recall Lemma 5.3 to estimate the expectation 
appearing in the statement of the lemma for large N. For small TV, it is 
enough to bound D+(t) by a Poisson point process. □ 

Lemma 5.6. For every T > and 5 > 0, 

rt+E 



lim lim supP^jv 



sup / £ N{U0) + Ul)}ds>S 

0<t<TJt 



0. 
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Proof. Recall the definition of the martingale M+(t) introduced in the 
beginning of the proof of Lemma 4.3. With this notation, in order to prove 
the lemma we need to show that 



lim limsupjr N 

£^0 iV^oo ^ 



sup {Mj v (t + e) 

0<t<T 



M^(t)}>5 



lim lim sup _ 



N 
N 



sup {D*(t + e)-D I *(t)}>6 



0<t<T 



0. 



0. 



and a similar statement with M^(t), D^(t) in place of M+ (*), £>+ (i). 

The martingale part is easy. By Doob's inequality and the explicit formula 
for the quadratic variation of M^(t) given at the beginning of the proof of 
Lemma 4.3, the probability which needs to be estimated is less than or equal 
to 



sup \M?(t)\>6/2 

0<t<T+e 



< W N 



T+e 



By the previous lemma, this expression vanishes as N j oo for any e > 0, 
5>0. 

On the other hand, the jump part has been estimated just after the proof 
of Lemma 4.5. □ 

We conclude this section with the proof of Proposition 4.2 which relies 
on the following lemma. For a subset / of R, denote by 1{/} the indicator 
function of /. 



Lemma 5.7. Fix T > and let [i N be a probability measure satisfying 
assumption (HI). There exist a finite constant Aq and a strictly positive 
continuous function Rt '■ (— oo, — Aq] U [Aq, oo) — ► R+ such that 



lim sup P^at 

TV^oo 



\(n?,l{[a,b]})\< I R T (u)du 







for all <t <T and all intervals [a, b] such that a > Aq or b < —Aq. 

Proof. Fix T > 0, < t < T and let \D N \(t) = D%(t) + D?(t) be the 
total number of particles which left the system before time t divided by N. 
By Lemma 5.4, there exists A\ > such that 

,/v r^^KT) > A x 



lim sup t „ 



0. 



Fix such Ai > and couple the £ process with a £ process as in Lemma 5.1 
with the additional property that (^-particles which jump to the interval 
{— A\ N, . . . , A\N} are removed. In particular, all ^-particles initially in this 
interval become instantaneously second-class particles. 
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On the set \D N \(T) < A\, we have that a particle is alive only if the 
particle Xj is alive and that the distance \Xj —Y^\, which does not depend 
on j for alive particles, is bounded by N\D N \(T). Therefore, if we fix an 
interval I = [a, b] with a > 3A 1 or b < -3A ± , on the set \D N \(T) < A 1 , 

(n N (Ct),l{I})=N- 1 £ Ct(x) = N- 1 J2nX{eNI} 1 

x/Nel j 

where the last summation is performed over all indices j corresponding to 
alive particles, and is bounded below by 

inf N- 1 ^2HY t ^N(v + I)}= inf (tt n (( t ) , l{v + I}) . 
\v\<Ai . M<Ai 

3 

In these formulas, v + I = {u + v :u 6 /} and NJ = {Nu :u£ J}. 

Since inf^,^^ (tt, 1{v+ J}) is a continuous function for the vague topology 
because all measures have density bounded by 1, by the hydrodynamic limit 
for the last expression converges in probability to 

inf (p t ,l{v + I})> / inf p(t,v)du, 

\v\<Aj J a \v—u\<A\ 

where p is the solution of the linear equation (5.2) with initial condition 
po and boundary condition p(t,±Ai) = 0. By the explicit formula (5.3) for 
the solution of (5.2), p(t,u) is smooth. Moreover, since po satisfies condi- 
tion (HI), p(t,u) is strictly positive. Therefore, for each u > 3Ai, 

R(T,u)= inf inf p(t,v)>0. 

0<t<T \v-u\<Ai 

This concludes the proof of the lemma. □ 

Proof of Proposition 4.2. In Lemma 5.7 we proved the proposition 
for intervals far from the origin. We estimate now the density on intervals 
close to the origin repeating the same argument presented in the proof of 
Lemma 5.7 and using the fact that the total number of particles which leave 
the system in a small time interval cannot be too large. 

Fix T > and recall from (4.6) that for each 5 > there exists e = e(5) > 
such that 



limsupP^v 



sup \D»\{t + e)-\D"\(t)>6 

0<t<T 



0. 



We shall estimate without loss of generality the density on M+. Fix a > 0. 
Let 5 = a/3, s = t — e{5) and denote by p N (s) the state of the £ process at 
time s. By Lemma 5.7, at time s, the particles' density is bounded below by 
a strictly positive function Rt: 

lim S u P /x JV ( S )[(7r JV ,l{/})<i? T (/)] = lim S upP^[{7rf <R T (I)} = 

N— >oo TV^oo 
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for all intervals / = [c, d] with c > Aq. Here, Rt(I) = Jj Rt( u ) du. 

Starting from [i N (s), we couple the £ process with a £ process as in the 
proof of Lemma 5.1 with the additional feature that (^-particles which reach 
the set {1, . . . , 5N} are killed. Following the argument presented in the proof 
of Lemma 5.7, on the set \D N \(t) - \D N \(s) < 5 we obtain that 

(7T N ^ t ),l{J})>^ N (Ct)Mv + J}), 

\v\<o 

for every interval J = [c, d] with c > a. The asymptotic behavior of the right- 
hand side of this inequality is given by the hydrodynamic limit of the £ 
process in the time interval [s,t]. The fact that we do not know the law of 
large number for the empirical measure at time s is not a problem. In fact, 
it is not difficult to show that tt n (Cr) is a tight sequence and that all limit 
points are concentrated on weak solutions of (5.2) with boundary condition 
p(r, S) = for s < r < t and on trajectories ir which at time s are bounded 
below by Rt on the interval f^O) 00 )- By monotonicity of weak solutions of 
(5.2), at time t the empirical measure has a density bounded below by the 
solution of (5.2) with initial condition Rt1{[Aq, oo)}. Therefore, in the limit 
N | oo, the right-hand side of the previous equation is bounded below by 

f d 

/ inf p(e,v)du, 

Jc \v—u\<8 

where p is the solution of (5.2) with initial condition pq = Rt1{[Aq, oo)} 
and boundary condition p(r,5) = for < r < e. Here again, the explicit 
formula (5.3) for the solution of (5.2) shows that the continuous function 
inf |„_ u i<j p(s, v) is strictly positive on [a,oo). This concludes the proof of 
the proposition. □ 
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